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Abstract. The 0/H metallicity distribution of different 
samples of planetary nebulae in the bulge of the Milky 
Way and M31 are compared. 0/H abundances are con- 
verted into [Fe/H] metallicity by the use of theoretical 
[O/Fe] X [Fe/H] relationships both for the bulge and the 
solar neighbourhood. It is found that these relationships 
imply an offset of [Fe/H] abundances by a factor up to 
0.5 dex for bulge nebulae. Systematic errors in the 0/H 
abundances as suggested by some recent recombination 
line work, ON cycling and statistical uncertainties are un- 
able to explain the observed offset, suggesting that the 
adopted relationship for the bulge probably overestimates 
the oxygen enhancement relative to iron. 
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1. Introduction 

Recent chemical evolution models usually predict different 
relationships between the [a-elements/Fe] and the metal- 
licity as measured by the [Fe/H] ratio for the different 
phases that comprise the Gala xy, na mely the disk, bulge 
and halo (see for example Pagel 1997 ). These relationships 
basically reflect the rate at which these elements are pro- 
duced in different scenarios, being usually faster in the 
bulge and halo compared to the disk in the framework of 
an inside-out model for Galaxy formation. 

Metallicities of bulge stars are poorly known compared 
with disk objects, as only limited samples of well measured 
stars are available. As a conclusion, the derived metallic- 
ity distribution and the corresponding ratios between the 
a-elements and metalUicity are not well known. In this 
work, a sample of bulge planetary nebulae (PN) with rel- 
atively accurate abundances is used to shed some light on 
the [O/Fe] x [Fe/H] relationship adopted for the bulge. It 
will be shown that this relation probably exagerates the 
amount of oxygen produced at a given metallicity. 
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Fig. 1. Lower panel: metallicity distribution of bulge PN 
from Cuisinier et al. (1999, thick line) and Ratag et al. 
(1997, thin line). Top panel: The same for the PN in the 
bulge of M31 (Jacoby & Ciardullo 1999) 



2. Metallicity distribution of bulge PN 

Many bulge, or type V PN (Maciel |l989| ) are known, but 
only recently have accurate abundances been obtained. 
Rece nt wo rk by Cuisinier et al. ( 1999 ) and Costa & Ma- 
ciel ( 1999 ) has led to He, O, N, Ar and S abundances for 
about 40 bulge PN, with an uncertainty comparable to 
disk objects, namely up to 0.2 dex for 0/H. The bulge 
0/H abundances are generally comparable with those of 
the disk, and the 0/H, Ar/H and S/H ratios can be higher 
than the disk counterparts even though very metal rich PN 
are missing in the bulge. Since underabundant nebulae are 
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also present, these results suggest that the bulge contains 
a mixed population, so that star formation in the bulge 
spans a wide time interval. 

Chemical abundances of PN in the bulge of M31 have 
been recently studied by Jacoby & CiarduUo ( |1999 ), who 
have included in their analysis some results by Stasihska 
et al. ( |1998D and Richer et al. ( |1999D . Fig. g shows a com- 
parison of this sample (42 PN, top panel) with the galac- 
tic bulge objects from Ratag et al. ( 1997|) (103 PN, lower 
panel, thin line) and Cuisinier et al. ( 1999| ) (30 PN, lower 
panel, thick line). The oxygen abundance is given in the 
usual notation, e(0) = log(0/H) -t- 12. It can be seen that 
all the 0/H abundance distributions are similar, peaking 
around 8.7 dex, and showing very few if any super metal 
rich objects with supersolar abundances. 

In order to compare the PN metallicity distribution 
with the stellar distributions, it is necessary to convert the 
measured nebular 0/H abundances into the usual [Fe/H] 
metallicities relative to the Sun. Direct measurements are 
of limited usefulness, for two main reasons. First, a very 
small number of planetary nebulae have measured iron 
lines, due to their weakness and the relatively large dis- 
tances of the nebulae, so that the derived values are more 
uncertain than the usual O, N or Ar abundances. Sec- 
ond, all available measurements indicate a strong deple- 
tion usually attributed to grain formation, so that the 
measured Fe abundances should be considered as lower 
limits to the total abundances at the times of formation 
of the PN progenitor stars. Both these aspects are illus- 
trated in Fig. m where we plot the [Fe/H] abundances 
against oxygen for a group of four disk planetary nebulae 
from Perinotto et al. ( 1999 , squares with e rror bars) and 
one object by Pottasch & Beintema ( 1999 , triangle with 
error bars). In all conversions we have used the solar iron 
abundance e(Fe)0 = 7.5 (Anders & Grevesse 1989 ) and 
average 0/H error bars. 

A better way to convert 0/H abundances into [Fe/H] 
metallicities is to use theoretical [0/Fe] x [Fe/H] relation- 
ships, such as those recently derived by Matteucci et al. 
( |l999|) both for the solar neighbourhood and the galactic 
bulge. The corresponding relations are more suitably plot- 
ted in Fig. I using e(0)o = 8.9 (Anders & Grevesse [l989[ ) 
both for the bulge (solid line) and the solar neighbourhood 
(dotted line). 

Taking into account the relations shown in Fig. H, the 



0/H distribution of bulge PN from Cuisinier et al. ( [1999D 
can be converted into a [Fe/H] distribution, as shown in 
Fig. y. Again the solid line corresponds to the [0/Fe] x 
[Fe/H] relationship for the bulge, and the dotted line was 
obtained using the solar neighbourhood relation shown in 
Fig. g. As a comparison, the dashed histogram in Fig. ^ 
shows the metallicity distribution of bulge K gia nt stars 
in Baade's Window by Mc William & Rich ( |l994| ). It can 
be seen that the PN metallicity distribution looks simi- 
lar to the K giant distribution if the solar neighbourhood 
[0/Fe] X [Fe/H] relation is adopted, but when the bulge 
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Fig. 3. Metallicity distribution of bulge PN using the 
[0/Fe] X [Fe/H] relations from JVIatteucci et al. (1999) 
for the bulge (solid line) and solar neighbourhood (dotted 
line). The dashed histogram shows the distribution of K 
giants from JVIcWiUiam & Rich (1994). 



relation is taken into account the derived distribution is 
displaced towards lower metallicities by roughly 0.5 dex. 
The PN samples were carefully selected not to include 
disk objects (see a detailed discussion in Cuisinier et al. 
1999), and their progenitor stars are clearly not massive 



enough to appreciably change their initial 0/H compo- 
sition. Therefore, the bulge PN metallicity distribution is 
expected to be similar to that of the K giants. This conclu- 
sion is strenghtened by the earlier measurements of Sadler 
et al. ( 1996 ) for JVIg and by the results recently presented 
by Feast (200C), who has shown that the metallicity distri- 
bution of bulge IVIira variables peaks around [Fe/H] ~ 0, 
which is again higher by roughly 0.5 dex than the bulge 
PN shown by the solid histogram in Fig. 0. Since both 
classes of objects are basically the offspring of the evolu- 
tion of intermediate to low mass stars, their metallicity 
distributions should also be similar. These objects may be 
located within a certain distance from the galactic centre, 
but, as discussed by Frogel (1999), any gradient between 
Baade's Window and the Centre must be small, amount- 
ing up to a few tenths of a dex. 

3. Discussion 

Several reasons can be considered in order to explain the 
discrepancy between the bulge PN distribution shown in 
Fig. I and the K giants of McWilliam & Rich ( [19941 ): 
(i) systematic errors in the 0/H abundances of planetary 
nebulae, leading to mild to strong underestimates of this 
quantity; (ii) ON cycling in the PN progenitor stars, which 
would lead to a depleted 0/H ratio; (iii) Statistical uncer- 
tainties due to the fact that the considered samples are 
relatively small and probably incomplete, and (iv) uncer- 



W.J. Maciel: Metallicity distribution of bulge PN 






o 




, 1 . 1 1 . 1 1 1 1 . 1 

+ 


+ 


+ 


1 1 1 


> 






Solar neighbourhood + . 


+ 




Bulge 


- 


^ 


- 


+ _,-•■ 


/ 








- 


' 


" 










/ k 










- 




■i. 


< \ 
















-■■;' 1 , , , /^\ 


- 








1 , , 




t^ 


1 




1 



7.5 8 8.5 

log (0/H) + 12 



9.5 



Fig. 2. [Fc/H] X log 0/H + 12 relation for the galactic bulge (solid line) and solar neighbourhood (dotted line) from 
the [0/Fe] x [Fe/H] relations of Matteucci et al. (1999). Also shown are: Four disk PN from Perinotto et al. (1999, 
squares with error bars), a disk PN from Pottasch & Beintema (1999, triangle with error bars), 6 bulge giants from 
McWilliam & Rich (1994, crosses), and a bulge star from Barbuy (1999, asterisk). 



tainties in the adopted [0/Fe] x [Fe/H] relationships. Let 
us briefly consider each of these possibilities. 

First, some recent work has raised the possibility that 
the 0/H abundances of planetary nebulae may be under- 
estimated, in view of the discrepancies between the forbid- 
den line and recombination line values for a number of ob- 
jects (Mathis & Liu 1999 , Liu & Danziger 1993| ). In prin- 



ciple, that would be applied to all PN samples con sidere d 
in this work, namely, the bulge PN of Ratag et al. (1997), 
Cuisinier et al. (1999) and the objects in the bulge of M31 
of Jacoby & CiarduUo ( 1999| ). There are no detailed pub- 
lished (under)abundance analyses for a large number of 
objects, but in order to be able to explain the discrep- 
ancy shown in Fig. ||, the 0/H abundances should have to 
be underestimated by about 0.5 dex, or a factor 3. This 
factor should also be applied to disk PN, as their abun- 
dances are derived using basically the same methods used 
for bulge PN. However, Cuisinier et al. (1999) have shown 
that their metallicity distribution for bulge PN is simi- 
lar to the correspo nding distribution of disk PN given by 
Maciel & Koppen (1994). This is confirmed by the recent 
(althoug h sma ller) sample of disk PN used by Maciel & 
Quireza ( 1999 ) to study radial abundance gradients in the 
galactic disk. This sample includes 128 disk PN, and an 
application of the solar neighb ourho od [0/Fe] x [Fe/H] re- 
lationship of Matteucci et al. ( 1999 ) (dotted line in Fig. ||) 
produces a [Fe/H] distribution peaked at [Fe/H] ~ —0.3. 
This is very similar to the recent G-dwarf metallicity dis- 
tribution of the solar neighbourhood by Rocha-Pinto & 



Maciel (1996) within the usually adopted uncertainties in 
the PN abundances of up to 0.2 dex for oxygen. Therefore, 
if we were to apply a correction factor of about 0.5 dex 
to the 0/H abundances, the derived disk PN distribution 
would imply a very large number of extremely metal-rich 
objects, whose nature it would be very difficult to explain. 
Moreover, the PN distribution would peak at about 9.2 dex, 
which means that most disk PN would be more oxygen rich 
than the Sun by about a factor 2. Even considering that the 
populations of PN and G-dwarfs have somewhat different 
ages, no known age-metallicity relation would be sufhcient 
to explain such large overabunda nce (se e for example the 
age metallicity relation by T warog 1980 or the recent work 
by Rocha-Pinto et al. 1999 and references therein). As a 
conclusion, any underestimate of the 0/H abundances in 
PN would have to be much lower than 0.5 dex, and could 
be well accommodated within an average uncertainty of 
0.2 dex. Of course, this does not exclude the possibility 
that some individual nebulae may have strongly under- 
estimated abundances, but the average factor is probably 
much lower than 0.5 dex. Further work is needed, possibly 
including a sizable sample of disk PN with well measured 
forbidden-line and recombination- line abundances. 

Second, the possibility of ON cycling has often been 
mentioned for planetary nebulae (see for example Maciel 
1992), basically due to some anticorrelation in the N/0 
ratio compared with He/H for disk PN. However, this 
phenomenon is only expected to occur in Type I PN (cf. 
Peimbert 1978), which are formed mainly from the higher 
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mass progenitor stars. These objects have an excess of 
He/H and/or N/0, and their oxygen abundances are in 
fact shghtly lower than t he mo st common Type II objects 
(see for example Maciel 2000 ). However, the amount by 
which the 0/H ratio is decreased is very small (roughly 
0.1 dex), and Type I objects are explici tly ex cluded from 
the disk sample of Maciel & Quireza ( 199£ ). Regarding 
the bulge PN, almost all objects show no trace of He/H or 
N/0 enrichment, so that this possibility cannot be used 
to explain the discrepancy shown in Fig. 0. 

Thirds statistical uncertainties are more difficult to an- 
alyze, since the considered samples are relatively small and 
may be affected by observational selection effects. How- 
ever, all distributions can be understood in terms of gen- 
eral models for chemical evolution of the Galaxy, and the 
similarities of the metallicity distributions of different ob- 
jects such as G-dwarfs and disk PN, or bulge PN, Mira 
variables and K giants fit rather nicely in the framework 
of galactic evolution, within the uncertainties of the de- 
rived abundances. Considering the bulge PN in particular, 
the similarity of the three different distributions shown in 
Fig. n^ is striking, even though they reflect different sam- 
ples using different techniques. Small differences s uch as 
the lack of oxygen rich PN in the Cuisinier et al. ( 1999 ) 
sample may be explained by the smaller size of this sam- 
ple. However, we are interested in very broad aspects of 
these distributions, and not in the detailed behaviour at 
a given metallicity of range of metallicities, so that it is 
unlikely that statistical uncertainties might produce the 
discrepancy shown in Fig. ||. 

Fourth, we are left with the possibility that the bulge 
[ 0/Fe ] X [Fe/H] relation as given by Matteucci et al. 
( 1999|) or, equivalently, the [Fe/H] x log 0/H relation 
given by the solid line in Fig. g might be responsible for all 
or most of the discrepancy in the metallicity distributions 
shown in Fig. ^. Such relation assumes a faster evolution 
during the bulge formation, so that at a given metallicity 
the [0/Fe] ratio is higher than in the solar neighbour- 
hood. Although this is correct in principle, the present 
results suggest that the amount of oxygen produced has 
been overestimated, leading to an excess of 0/H for a 
given metallicity. In fact, the earlier models of Matteucci 
& Brocato ( |1990D predict a lower [0/Fe] enhancement, 
producing a better agreement with the present results. On 
the other hand, there are some independent evidences that 
the theoretical [0/Fe] in the bulge may be overestimated. 
A recent analysis by Barbuy and collaborators for one star 
in the b ulge globular cluster NGC 6528 (see for example 
Barbuy |l999|) shows that [0/Fe] ~ 0.35 for a metallicity 
[Fe/H] ~ -0.6, so that [0/H] ~ -0.25 and log 0/H + 12 
~ 8.65. This object is also shown in Fig. g (asterisk) and 
is closer to the dotted curve than to the bulge relation 
(solid line). Also, for all six bulge giants with measur- 
able [01] features in the Mc William & Rich (|1994[ ) sample 
the [0/Fe] ratio is lower by 0.2 — 0.7 dex than predicted 
by the [0/Fe] x [Fe/H] relationship of Matteucci et al. 



( 1999D . These stars are also shown in Fig. || (crosses) , and 
are clearly located to the left of the bulge curve, even al- 
lowing for some uncertainty in the 0/H abundances. The 
average difference is about 0.5 dex, which is just what is 
needed to eliminate the offset of the bulge PN metallicity 
distribution shown in Fig. |[ Therefore, it can be con- 
cluded that the theoretical [0/Fe] x [Fe/H] relation for 
the bulge probably overestimates the oxygen enhancement 
relative to iron by 0.3 to 0.5 dex, at least for metallicities 
[Fe/H] > -1.5 dex. 
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